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These days, as high energy particle colliders become unavailable for testing speculative 
theoretical ideas, physicists are looking to other environments that may provide extreme 
conditions where theory confronts physical reality. One such circumstance may arise at high 
temperature T, which perhaps can be attained in heavy ion collisions or in astrophysical 
settings. It is natural therefore to examine the high-temperature behavior of the standard 
model, and here I shall report on recent progress in constructing the high-T limit of QCD. 

My presentation will be unified by the theme of screening, a familiar phenomenon 
in electro dynamical plasmas. I shall explore how similar effects can be described in QCD 
at a sufficiently high temperature (above the putative confinement - deconfinement phase 
transition) so that we may speak of unconfirmed quarks and gluons forming a plasma. But 
first let me review briefly the screening phenomena in plasmas of electromagnetically charged 
particles. We begin with Poisson's equation, which relates the scalar electric potential to 
a charge density p. 

- w = p 

For the charge density we take a statistical distribution of positive-charged (+q) and negative- 
charged (—q) particles, each carrying the energy ±g0, respectively, and described by the same 
density n. Then 

p = n (qe~ q(t>/T - qe q(t>/T ) 
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For large T, this becomes 



P 



(large T) 



2nq 2 <p/T 



so that the Poisson equation reads 



-V 2 + 



2nq 2 



T 



) 



= 



Evidently, a screening mass oc ( j 2 has been induced for the electric potential 0; the inverse 
is called the Debye screening length. Again at high T and for a relativistic plasma, one 
expects n ~ (1/volume) ~ T 3 , hence the induced electric screening mass is 



We shall see a similar result emerging in the non-Abelian theory as well. Note that Debye 
screening occurs for the electric (temporal) component of the gauge potential. There is no 
electro dynamical magnetic screening, because there are no magnetic resources. 



In the non-Abelian theory, the corresponding equation involves the covariant divergence. 



(Here g is the gauge coupling constant.) So the issue of magnetic sources is not so clear 
in the Yang-Mills case, and one of the topics that we shall address later is whether in the 
non-Abelian theory there exists magnetic screening. 

The above argument - it is essentially Debye's - makes little use of field theoretical 
formalism. But to carry through analogous calculations in the standard model, we shall 
begin with quantum field theory. Let me explore how finite-temperature calculations are 
performed in that context. 

When studying a field theory at finite temperature, the simplest approach is the so- 
called imaginary-time formalism. We continue time to the imaginary interval [0, 1/iT] and 
consider bosonic (fermionic) fields to be periodic (anti-periodic) on that interval. Perturba- 
tive calculations are performed by the usual Feynman rules as at zero temperature, except 
that in the conjugate energy-momentum, Fourier-transformed space, the energy variable 
p° (conjugate to the periodic time variable) becomes discrete - it is 2imT, (n integer) for 
bosons. From this one immediately sees that at high temperature - in the limiting case, at 
infinite temperature - the time direction disappears, because the temporal interval shrinks 
to zero. Only zero-energy processes survive, since "non-vanishing energy" necessarily means 
high energy owing to the discreteness of the energy variable p° ~ 2irnT, and therefore all 
modes with n ^ decouple at large T. In this way a Euclidean three-dimensional field 
theory becomes effective at high temperatures and describes essentially static processesJil 

Let me repeat in greater detail. Finite-T, imaginary-time perturbation theory makes 
use of conventional diagrammatic analysis in "momentum" space, with modified "energy" 
variables, as indicated above. Specifically a spinless boson propagator is 



m oc \q\T 



VB = O 



V • B a = gf abc A b ■ B c 
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while a spin-| fermion propagator reads 



S(p) 



Po = m(2n + 1)T 



7 Po — 1 • P — m 



The zero-temperature integration measure / becomes replaced by 

iTJ2™=-oc I J0js- Thus it is seen that Bose exchange between two 0(g) vertices contributes 
iT E£L-oo J (fjr 9 - 47r 2 n 2 T 2_pii_ m 2 9 where g is the coupling strength. In the large T limit, 
all n 7^ terms (formally) vanish as T~ l and only the n = term survives. One is left 
with f -£^Lg^/T .2} 2 gVT. This is a Bose exchange graph in a Euclidian 3-dimensional 

theory, with effective coupling gyT. Similar reasoning leads to the conclusion that fermions 
decouple at large T. 

While all this is quick and simple, it may be physically inadequate. First of all, fre- 
quently one is interested in non-static processes in real time, so complicated analytic con- 
tinuation from imaginary time needs to be made before passing to the high-T limit, which 
in imaginary time describes only static processes. Also one may wish to study amplitudes 
where the real external energy is neither large nor zero, even though virtual internal energies 
are high. 

Another reason that the above may be inadequate emerges when we consider massless 
fields (such as those that occur in QCD). We have seen that the n = mode leaves a 
propagator that behaves as ^ when mass vanishes, and a phase space of d 3 p. It is well known 
that this kind of kinematics at low momenta leads to infrared divergences in perturbation 
theory even for off-mass-shell amplitudes - - Green's functions Jn massless Bosonic field 
theories possess infrared divergences in naive perturbation theory.H Since physical QCD does 
not suffer from off-mass-shell infrared divergences, perturbation theory must be resummed. 

A final shortcoming of the above limiting procedure is that it is formal: the limit is 
taken before the integration/summation is carried out. But the latter need not converge 
uniformly; indeed owing to ultraviolet divergences, it may not converge at all and must be 
renormalized. As a consequence the n ^ contributions in single Boson exchange graphs 
may not decrease as T _1 . 

Thus the formal arguments for the emergence of a 3-dimensional theory at high-T 
need be re-examined for QCD. Nevertheless, even if unreliable, the arguments alert us 
to the possibility that 3-dimensional field theoretic structures may emerge in the high-T 
regime. Indeed this occurs, although not in a direct, straightforward fashion; this will be 
demonstrated presently. 

Here is a graphical argument to the same end discussed above: viz. The need to resum 
perturbation theory. Consider a one-loop amplitude Ui(p), 




given by the graph in the figure. 
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= J dk Ii(p, k) 

Compare this to a two- loop amplitude n 2 (p), 

n 2 (p) = J dkl 2 {p,k) , 
in which 111 is an insertion, as in the figure below. 



n 2 (p) = 



= J dk h(p,k) 

Following Pisarski,i we estimate the relative importance of H2 to 111 by the ratio of their 
integrands, 

112 12 2 n x (fc) 
n x ~ h 9 k 2 ' 

Here g is the coupling constant, and the k 2 in the denominator reflects the fact that we are 
considering a massless field, as in QCD. Clearly the k 2 — > limit is relevant to the question 
whether the higher order graph can be neglected relative to the lower order one. Because 
one finds that for small k and large T, Ili(fc) behaves as T 2 , the ratio H2/H1 is g 2 T 2 /k 2 . As 
a result when k is O(gT) or smaller the two- loop amplitude is not negligible compared to 
the one-loop amplitude. Thus graphs with "soft" external momenta [O(gT) or smaller] have 
to be included as insertions in higher order calculations. 

A terminology has arisen: graphs with generic/soft external moment [O(gT) where g 
is small and T is large] and large internal momenta [the internal momenta are integration 
variables in an amplitude; when T is large they are 0(T), hence also large] are called "hard 
thermal loops. "IS Much study has been expended on them and finally a general picture has 
emerged. Before presenting general results, let us look at a specific example — a 2-point 
Green's function. 

It needs to be appreciated that in the imaginary-time formalism the correlation func- 
tions are unique and definite. But passage to real time, requires continuing from the integer- 
valued "energy" to a continuous variable, and this cannot be performed uniquely. This 
reflects the fact that in real time there exists a variety of correlation functions: time ordered 
products, retarded commutators, advanced commutators, etc. Essentially one is seeing the 
consequence of the fact that a Euclidean Laplacian possesses a unique inverse, whereas giv- 
ing an inverse for the Minkowskian d'Alembertian requires specifying temporal boundary 
conditions, and a variety of answers can be gotten with a variety of boundary conditions. 
Thus, when presenting results one needs to specify precisely what one is computing. 

We shall consider a correlation function for two fermionic currents, in the 1-loop ap- 
proximation. 

W(x,y) = -i (nx)f(y)) 
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d A k 



-ik(x-y) 



W u (k) 



The QCD result differs from the QED result by a group theoretical multiplicative factor, so 
we present high-T results only for the latter, in real-time, and consider the time-ordained 
product 11%? as well as the retarded commutator H 1 ^ ■ 

IP" possesses a real and an imaginary part. It is found that at large T, the real parts 
of lijf and 11^ coincide. 



T 2 T 2 k 2 
-ReIP"(£;) = —P% v + 



\k l 



1 + 4^ In 
2 k 



k° - \k\ 



k° + \k\ 



where the projection operators are 



if \x or v = 

5 lj — k' l k j /\k 2 \ otherwise 



For the imaginary part, which is present only for space-like arguments, different expressions 
are found. 



-Imn^(ifc) = -ixp^ik) 



Ifcl 3 2 



T 2 9(-k 2 



-Imn£?(Jfe) 



nk 2 



k° 
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9(-k 2 ) 



3 1 + 2 2 . 



A unified presentation of these formulas is achieved in a dispersive representation. For 
the retarded function this reads 



nST(fc) = n& fl (*) + 

while the time-ordered expression is 

nT(*) = n& B (*) + 



/ P^(fco,fe) 



2ni 



p^(k ,k) 



v °k' -k -ie ' e k °/ T - , 
The dispersive expressions may also be used to give the imaginary-time formula. 



imaginary 
time 



(*) 



dk 



k' — 2mnT 



In all the above formulas, ^sub * s a rea l subtraction term. 

Note that a universal statement about high-T behavior can be made only for the 
absorptive part : it is 0{T 2 ). This also characterizes 11^, but lijf possesses an additional 
0(T 3 ) imaginary part, which is seen to arise from the additional term in Hjf involving the 
bosonic distribution function efco / T _ 1 - Finally, the II^ a g inary amplitude has a temperature 
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behavior determined by its external "energy" = 2irinT. If this is replaced by a fixed ko 
(T-independent) or if only the n = mode is considered, then one may assign an 0(T 2 ) 
behavior to this quantity as well. 

In conclusion, we assert that the 2-point correlation function behaves as 0(T 2 ), where 
it is understood that this statement is to be applied to the retarded amplitude, ox to the 
imaginary time amplitude with its "energy" argument continued away from 2-KinTB 

Similar analysis has been performed on the higher-point functions and this work has 
culminated with the discovery (Braaten, Pisarski, Frenkel, Taylor)BoQ of a remarkable sim- 
plicity in their structure. To describe this simplicity, we do not discuss the individual n-point 
functions, but rather their sum multiplied by powers of the vector potential, viz. we consider 
the generating functional for single-particle irreducible Green's functions with gauge field ex- 
ternal lines in the hard thermal limit. (Effectively, we are dealing with continued imaginary- 
time amplitudes.) We call this quantity r H TL(^4) and it is computed in an SU(N) gauge 
theory containing Np fermion species of the fundamental representation. Thtl is found (i) 
to be proportional to (N + ^Np), (ii) to behave as T 2 at high temperature, and (iii) to be 
gauge invariant. 



A A 

A„ A 2, <n > 



A~ ' A 



i 
A 



(N+lN F ) 9 -^T HTh (A) 
IZn 



Tn^iU- 1 AU + U- 1 dU) = r HTL (A) 

(Henceforth g, the coupling constant, is scaled to unity.) A further kinematical simplification 
in Thtl has also been established. To explain this we define two light-like four-vectors Q± 
depending on a unit three-vector q, pointing in an arbitrary direction. 



q ■ q = 1 , Q'±Q ±M = , Q^Q TI1 = 1 
Coordinates and potentials are projected onto Q±. 

X ± EE X,Q^ , 8± EE Ql^jJ > A ± EE 

The additional fact that is now known is that (iv) after separating an ultralocal contribution 
from Thtl, the remainder may be written as an average over the angles of q of a functional 
W that depends only on A + ; also this functional is non-local only on the two-dimensional 
x plane, and is ultralocal in the remaining directions, perpendicular to the x plane. ["Ul- 
tralocal" means that any potentially non-local kernel k(x,y) is in fact a ^-function of the 
difference k(x,y) oc 8(x — y).] 

T HTL (A) = 2vr | d 4 x A a (x)A a (x) + J dn 4 W(A + ) 
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These results are established in perturbation theory, and a perturbative expansion of W(A + ), 
i.e. a power series in A + , exhibits the above mentioned properties. A natural question is 
whether one can sum the series to obtain an expression for W(A + ). n 
Important progress on this problem was made when it was observed (Taylor, Wongp 
that the gauge-invariance condition can be imposed infinitesimally, whereupon it leads to a 
functional differential equation for which is best presented as 



d 



dx+ 5AI 



W{A + ) + \ / <fx A\(x)A b Ax) 







—dx~ 



abc <\b 



+ f aoc A 
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5A C , 



W{A+) + \ / d 4 x A d Jx)A d Jx) 



In other words we seek a quantity, call it 

S(A+) = W(A+) + \J *xA%(x)A%(x) , 

which is a functional on a two-dimensional manifold {x + , x~}, depends on a single functional 
variable A + , and satisfies 



<9i— S - d 2 A? + f abc A b —S 



c -S -d 2 A°,+ f abc A\- 
8AI 1 J 1 8A C 1 







x 



'2" = -af 



A\ = A\ 



Another suggestive version of the above is gotten by defining A\ = j^. 



d x A a 2 - d 2 Al + f abc A\A c 2 = 

To solve the functional equation and produce an expression for W(A+), we now turn to a 
completely different corner of physics, and that is Chern-Simons theory at zero temperature. 

The Chern-Simons term is a peculiar gauge theoretic topological structure that can be 
constructed in odd dimensions, and here we consider it in 3-dimensional space-time. 



J cs « J d 3 x e a ^ Tr {d a A p A 1 + § A a ApA, 



This object was introduced into physics over a decade ago, and since that time it has been put 
to various physical and mathematical uses. Indeed one of our originally stated motivations 
for studying the Chern-Simons term was its possible relevance to high-temperature gauge 
theory.Q Here following Efraty and Nairja we shall employ the Chern-Simons term for a 
determination of the hard thermal loop generating functional, Thtl- 

Since it is the space-time integral of a density, Jcs m ay be viewed as the action for a 
quantum field theory in (2+l)-dimensional space-time, and the corresponding Lagrangian 
would then be given by a two-dimensional, spatial integral of a Lagrange density. 



J cs oc J dt L cs 

L cs oc J d 2 x (A a 2 A1 + A a F? 2 ) 
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I have separated the temporal index (0) from the two spatial ones (1,2) and have indicated 
time differentiation by an over dot. F£ 2 is the non-Abelian field strength, defined on a 
two-dimensional plane. 

= d x A a 2 - d 2 A a 1 + f ahc A\A\ 
Examining the Lagrangian, we see that it has the form 

L ~ pq — A H (p, q) 

where A 2 plays the role of p, A\ that of q, F® 2 is like a Hamiltonian and Aq acts like the 
Lagrange multiplier A, which forces the Hamiltonian to vanish; here Aq enforces the vanishing 

F° = 

The analogy instructs us how the Chern-Simons theory should be quantized. 

We postulate equal-time computation relations, like those between p and q. 



A°(r), A b 2 (r')\ = i5 ab 5(r-r') 

In order to satisfy the condition enforced by the Lagrange multiplier, we demand that F± 2 , 
operating on "allowed" states, annihilate them. 







This equation can be explicitly presented in a Schrodinger-like representation for the 
Chern-Simons quantum field theory, where the state is a functional of Af. The action of the 
operators A\ and A 2 is by multiplication and functional differentiation, respectively. 



A\ 



A a 2 



) 



) 



Al^f(A^) 
iSAt v 1J 



This, of course, is just the field theoretic analog of the quantum mechanical situation where 
states are functions of q, the q operator acts by multiplication, and the p operator by dif- 
ferentiation. In the Schrodinger representation, the condition that states be annihilated by 



pa 



'd^Al - d 2 A\ + j ahc A\A\ 
leads to a functional differential equation. 







1 5 



i 5A a 



d 2 A a 1 + f abc A b 1 



1 5 



i 8A\ 



= 



If we define S by ^ = e lS we get equivalently 
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This equation comprises the entire content of Chern-Simons quantum field theory. S is the 
Chern-Simons eikonal, which gives the exact wave functional owing to the simple dynamics 
of the theory. Also the above eikonal equation is recognized to be precisely the equation for 
the hard thermal loop generating functional, given above. 

Let me elaborate on the connection with eikonal- WKB ideas. Let us recall that in 
particle quantum mechanics, when the wave function ip{q) is written in eikonal form 

tj}{ q ) = e iS{q) 

then the WKB approximation to S(q) is given by the integral of the canonical 1-form pdq 

S(q) = / p(q')dq 



where p(q), the momentum, is taken to be function of the coordinate q, by virtue of satisfying 
the equation of motion. 



p(q) = pE - 2V(q) 

Analogously, in the present field theory application, the eikonal S(Ai) may be written as a 
functional integral, 

S{A{) = J 1 Al{A' x )VA'? 
where Af(Ai) is functional of A\ determined by the equation of motion 

d x A% - d 2 A\ + j ahc A\A\ = 

Since, by construction jjjp = A?,, it is clear that as a consequence S satisfies the required 

equation. However, we reiterate that in the Chern-Simons case there is no WKB approxi- 
mation: everything is exact owing to the simplicity of Chern-Simons dynamics. 

The gained advantage for thermal physics is that "acceptable" Chern-Simons states, 
i.e. solutions to the above functional equations, were constructed long ago,H and one can 
now take over those results to the hard thermal loop problem. One knows from the Chern- 
Simons work that \1/ and S are given by a 2-dimensional fermionic determinant, i.e. by the 
Polyakov-Wiegman expression. While these are not described by very explicit formulas, 
many properties are understood, and the hope is that one can use these properties to obtain 
further information about high-temperature QCD processes. We give two applications. 

The Chern-Simons information allows presenting the hard-thermal loop generating 
functional as 

r H TL = \jdn 4 [A a + A a _ + S(A+)+S{A.)] . 

Using the known properties of S, one can give a very explicit series expansion for Thtl in 
terms of powers of A 

1 f-m . . 1 



r HTL = jj T^AA + -J T^AAA + 
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where the non-local kernels are known explicitly. This power series may be used to 

systematize the resummation procedure for the pertubative theory. Here is what one does: 
perturbation theory for Green's functions may be organized with the help of a functional 
integral, where the integrand contains (among other factors) e* /QCD ^ A - ) where Iqcd is the 
QCD action. We now rewrite that as 



i 

e 



{/QCD(A)+^r H TL(A)-^r H TL(^)} 



where m = TJ N+ ^ F ' 2 _ Obviously nothing has changed, because we have merely added and 
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subtracted the hard-thermal-loop generating functional. Next we introduce a loop counting 
parameter I: in an /-expansion, different powers of I correspond to different numbers of loops, 
but at the end / is set to unity. The resummed action is then taken to be 



g^resumracd __ g 



One readily verifies that an expansion in powers of / describes the resummed perturbation 
theory, and this then represents the first application of the present Chern-Simons formalism. 

For a second application, we note that even though the closed form for Thtl is not very 
explicit, a much more explicit formula can be gotten for its functional derivative ^"a L . This 

may be identified with an induced current, which is then used as a source in the Yang- Mills 
equation. Thereby one obtains a non-Abelian generalization of the Kubo equation, which 
governs the response of the hot quark gluon plasma to external disturbances.^ 



2 

T77 

n f^ v = j u 



induced 



From the known properties of the fermionic determinant — hard thermal loop generating 
functional — one can show that j^ duced is given by 

Educed = / ^ [<T + («-(*) - M*)) + Q- («+(*) " } 

where a± are solutions to the equations 

<9 + a_ - d_A + + [A + , a_] = 
d + A_ - <9_a + + [a+, AJ\ = 

Evidently ji nduced , as determined by the above equations, is a non-local and non-linear func- 
tional of the vector potential A^. 

Ther^now have appeared several alternative derivations of the Kubo equation. Blaizot 
and IancutO have analyzed the Schwinger- Dyson equations in the hard thermal regime; they 
truncated them at the 1-loop level, made further kinematical approximations that are justi- 
fied in the hard thermal limit, and they too arrived at the Kubo equation. Equivalently the 
argument may be presented succinctly in the language of the composite effective action,0 
which is truncated at the 1-loop (semi-classical) level — two-particle irreducible graphs are 
omitted. The stationarity condition on the 1-loop action is the gauge invariance constraint 
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on Thtl- Finally, there is one more, entirely different derivation — which perhaps is the 
most interesting because it relies on classical physics.tj We shall give the argument presently, 
but first we discuss solutions for the Kubo equation. 

To solve the Kubo equation, one must determine a± for arbitrary A±, thereby obtaining 
an expression for the induced current, as a functional of A±. Since the functional is non- 
local and non-linear, it does not appear possible to construct it explicitly in all generality. 
However, special cases can be readily handled. 

In the Abelian case, everything commutes and linearizes. One can determine a± in 
terms of A±. 

^ A 

d T 

(Incidentally, this formula exemplifies the kinematical simplicity, mentioned above, of hard 
thermal loops: the nonlocality of \/d± lies entirely in the {x + ,x~} plane.) With the above 
form for a± inserted into the Kubo equation, the solution can be constructed explicitly. 
It coincides withthe results obtained by Silin long ago, on the basis of the Boltzmann- 
Vlasov equation.tj One sees that the present theory is the non-Abelian generalization of 
that physics; in particular m, given above, is recognized as the Debye screening length, 
which remains gauge invariant in the non-Abelian context. 

It is especially interesting to emphasize that Silin did not use quantum field theory in 
his derivation; rather he employed classical transport theory. Nevertheless, his final result 
coincides with what here has been developed from a quantal framework. This raises the 
possibility that the non-Abelian Kubo equation can also be derived classically, and indeed 
such a derivation has been given, as mentioned above. 

We now pause in our discussion of solutions to the non-Abelian Kubo equation in order 
to describe its classical derivation. 

Transport theory is formulated in terms of a single-particle distribution function / 
on phase space. In the Abelian case, / depends on position {x^} and momentum {p^} 
of the particle. For the non-Abelian theory it is necessary to take into account the fact 
that the particle's non-Abelian charge {Q a } also is a dynamical variable: Q a satisfies an 
evolution equation (see below) and is an element of phase space. Therefore, the non-Abelian 
distribution function depends on {x^}, {p^} and {Q a }, and in the collisionless approximation 
obeys the transport equation j-f = 0, i.e. 

df dx» df dp^ df dQ a _ Q 



dx^- dr dp^ dr dQ a dr 

The derivatives of the phase-space variables are given by the Wong equations, for a particle 
with mass \i. 

dx** p^ 
dr \i 

dr a dr 

_ jabc qc 

dr dr M 
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In order to close the system we need an equation for F^ v '. In a microscopic description (with 
a single particle) one would have {D^F^ p ) a = J dTQ a {r)^-^-8 A {x — x(r)^j and consistency 
would require covariant conservation of the current; this is ensured provided Q a satisfies the 
equation given above. In our macroscopic, statistical derivation, the current is given in terms 
of the distribution function, so the system of equations closes with 



(ZV^T = \d v dQQ a ff{x^Q) 



(One verifies that the current - the right side of the above - is covariatly conserved.) The col- 
lisionless transport equation, with the equations of motion inserted, is called the Boltzmann 
equation. The closed system formed by the latter supplemented with the Yang-Mills equa- 
tion is known as the non-Abelian Vlasov equations. To make progress, this highly non-linear 
set of equations is approximated by expanding around the equilibrium form for /, 

boson V 
fermion 

This^cpmprises the Vlasov approximation, and readily leads to the non-Abelian Kubo equa- 
tion.!^ 

One may say that the non-Abelian theory is the minimal elaboration of the Abelian 
case needed to preserve non-Abelian gauge invariance. The fact that classical reasoning 
can reproduce quantal results is presumably related to the fact that the quantum theory 
makes use of the (resummed) 1-loop approximation, which is frequently recognized as an 
essentially classical effect. Evidently, the quantum fluctuations included in the hard thermal 
loops coincide with thermal fluctuations. 

Returning now to our summary of the solutions to the non-Abelian Kubo equation 
that have been obtained thus far, we mention first that the static problem may be solved 
completely .EJ When the Ansatz is made that the vector potential is time independent, A± = 
A±(r), one may solve for a± to find a± = —A± and the induced current is explicitly computed 
as 

m 2 -jU / — m 2 v4°\ 

2 ^induced I g I 

This exhibits gauge-invariant electric screening with Debye mass m. One may also search for 
localized static solutions to the Kubo equation, but one finds only infinite energy solutions, 
carrying a point-magnetic monopole singularity. Thus there are no plasma solitons in high-T 
QCD.O Specifically, upon selecting the radially symmetric solution that decreases at large 
distances, there arises a magnetic monopole-like singularity at the origin. 

Much less is known concerning time-dependent solutions. Blaizot and IancuE-3 have 
made the Ansatz that the vector potentials depend only on the combination x ■ k, where 
k is an arbitrary 4- vector: A± = A±(x ■ k). Once again a± can be determined; one finds 

a± = ^^A T , and the induced current is computable. For k = Qj, where there is no space 

dependence (only a dependence on time is present) one finds 

^, = ( \ 
2 Jinduccd \-lm 2 AJ 
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More complicated expressions hold with general k. The Kubo equation can be solved nu- 
merically; the resulting profile is a non-Abelian generalization of a plasma plane wave. 

The physics of all these solutions, as well as of other, still undiscovered ones, remains 
to be elucidated, and I invite any of you to join in this interesting task. 

We see that Debye electric screening is reproduced in essentially the same form as 
in an Abelian plasma (to leading order). How about magnetic screening? It is important 
to appreciate that the above time- independent, space-independent induced current, with j 
proportional to A, does not describe magnetic screening because screening is determined 
by static configurations. Thus we conclude that the hard-thermal-loop limit of hot QCD 
does not show magnetic screening. Indeed it appears that if one proceeds perturbatively, 
beyond the resummed perturbation expansion of hard thermal loops, no direct evidence for 
magnetic screening can be found. 

However, there is indirect evidence: although the hard thermal loop resummation cures 
some of the perturbative infrared divergences, as one calculates to higher perturbative or- 
ders, they reappear essentially due to the non-linear interactions between electric (temporal) 
and unscreened magnetic (spatial) degrees of freedom as well as among the magnetic de- 
grees of freedom due to their self-interaction. (Such interactions are absent in an Abelian 
theory.) Consequently it is believed that non-perturbative magnetic screening arises in the 
non-Abelian theory, and it is recalled that, as mentioned in the Introduction, there is some- 
thing akin to a magnetic source in Yang-Mills theory. 

Another qualitative argument can be offered to make plausible the idea that a magnetic 
mass should arise. Although I have argued that high-temperature dimensional reduction from 
four to three dimensions can not be carried out reliably for a gauge theory, one may speculate 
that there is some truth in the idea, when restricted to magnetic (spatial) components of the 
non-Abelian potential. So one is led, as preliminary to studying the full QCD problem, to an 
analysis of three-dimensional Euclidean Yang-Mills theory at zero temperature. One quickly 
discovers that infrared divergences are present in perturbation theory for this model as well, 
so here again arises the question of a dynamically induced mass. In three dimensions, the 
coupling constant squared carries dimensions of mass. (Recall that in a high-temperature 

reduction is related to the four-dimensional coupling g by g^ = gVT.) Therefore 
it is plausible that three-dimensional Yang-Mills theory generates an 0(g? 3 \) mass, which 
eliminates its perturbative infrared divergences and suggests the occurrence of an 0(g 2 T) 
magnetic mass in the four- dimensional theory at high T. Unfortunately thus far no analysis 
of the three-dimensional Yang-Mills model has led to a proof of such mass generation. 

Since the mass is not seen in perturbative expansions, even resummed ones, one at- 
tempts a non-perturbative calculation, based on a gap equation. Of course an exact treat- 
ment is impossible; one must be satisfied with an approximate gap equation, which effectively 
sums a large, but still incomplete set of graphs. At the same time, gauge invariance should 
be maintained; gauge non-invariant approximations are not persuasive. 

Deriving an approximate, but gauge invariant gap equation is most efficiently carried 
out in a functional integral formulation. We begin by reviewing how a one-loop gap equation 
is gotten from the functional integral, first for a non-gauge theory of a scalar field tp, then 
we indicate how to extend the procedure when gauge invariance is to be maintained for a 
gauge field A^. 
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Consider a self-interacting scalar field theory (in the Euclidean formulation) whose 
potential V(ip) has no quadratic term, so in direct perturbation theory one may encounter 
infrared divergences, and one enquires whether a mass is generated, which would cure them. 

C = \d^d^ + V{<p) 
V(<p) = A 3 y? 3 + A 4 y? 4 + • • • 

The functional integral for the Euclidean theory involves the negative exponential of the 
action I — J C. Separating the quadratic, kinetic part of /, and expanding the exponential 
of the remainder in powers of the field yields the usual loop expansion. As mentioned earlier, 
the loop expansion may be systematized by introducing a loop-counting parameter £ and 
considering e~ j 11 ^^) •. the power series in £ is the loop expansion. To obtain a gap equation 
for a possible mass /x, we proceed by adding and subtracting 1^ = ^ J ip 2 , which of course 
changes nothing. 

Next the loop expansion is reorganized by expanding I + 1^ in the usual way, but taking 
—1^ as contributing at one loop higher. This is systematized as in the hard-thermal-loop 
application with an effective action, 1^ containing the loop counting parameter £, which 
organizes the loop expansion in the indicated manner: 

h = j (i(Ve<p) + i^Viip)) - UVIcp) 

An expansion in powers of £ corresponds to a resummed series; keeping all terms and setting 
I to unity returns us to the original theory (assuming that rearranging the series does no 
harm); approximations consist of keeping a finite number of terms: the 0(£) term involves 
a single loop. 

The gap equation is gotten by considering the self energy £ of the complete propagator. 
To one-loop order, the contributing graphs are depicted in the Figure. 




Self energy resummed to one-loop order. 

Regardless of the form of the exact potential, only the three- and four- point vertices 
are needed at one-loop order; the "bare" propagator is massive thanks to the addition of the 
mass term jl^Vltp) = UV; the last — /i 2 in the Figure comes from the subtraction of the 

same mass term, but at one-loop order: —I^(^/£ip) = —£\ J y? 2 '■ 

The gap equation emerges when it is demanded that E does not shift the mass ji. In 
momentum space, we require 
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E(p) =0 
\p 2 =— p? 




Graphical depiction of above equation. 

While these ideas can be applied to a gauge theory, it is necessary to elaborate them so 
that gauge invariance is preserved. We shall discuss solely the three-dimensional non-Abelian 
Yang-Mills model (in Euclidean formulation) as an interesting theory in its own right, and 
also as a key to the behavior of spatial variables in the physical, four- dimensional model at 
high temperature. 

The starting action I is the usual one for a gauge field. 

I = J d 3 x tr ^F i F i 

pi _ l Jjkp 
r ~ 2 l k 

While one may still add and subtract a mass-generating term J M , it is necessary to preserve 
gauge invariance. Thus we seek a gauge invariant functional of I^(A), whose quadratic 
portion gives rise to a mass. Evidently 

IM) = ~ y / d 3 x tr At Us -^).t ; • ... 

The transverse structure in the above equation guarantees invariance against Abelian gauge 
transformations; the question then remains how the quadratic term is to be completed in 
order that Ip(A) be invariant against non-Abelian gauge transformations. [In fact for the 
one-loop gap equation only terms through 0(A 4 ) are needed. 1 

A very interesting proposal for I^{A) was given by Nairt3El who also put forward the 
scheme for determining the magnetic mass, which we have been describing. By modifying in 
various ways the hard thermal loop generating functional (which gives a four-dimensional, 
gauge invariant but Lorentz non- invariant effective action with a transverse quadratic term), 
he arrived at a gauge and rotation invariant three-dimensional structure, which can be em- 
ployed in the derivation of a gap equation.^ 

Let me describe Nair's modification. Recall that the hard-thermal-loop generating 
functional, which I record here again, 

Thtl = \J dn 4 [Q^A«Al + S(Q»A,) + S(Q»A,)} 

* A gap equation for the full gauge-field propagator, rather than just for the mass, has been put forward and 
analyzed by Cornwall el al; see Phys. Lett. B153, 173 (1985) and Phys. Rev. D 34, 585 (1986). 
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is gauge invariant because Q± = (0, ±q) is light-like and S is the Chern-Simons eikonal. 
Before averaging over Qq one is dealirjffwith a functional of only A±\ after averaging all four 
components of enter. With Nair ,113119 we observe that another choice for Q± can be made, 
where those vectors remain light-like, but have vanishing time component. This is achieved 
when the spatial components of Q± are complex and of zero length; for example: 

Qt = (0,q) e = (0,q*) 

q = (— cos 9 cos ip — i sin ip, — cos 9 sirup + i cos ip, sin 9) = 9 + itp 

Evidently q 2 = 0, Q\ = and Q- = Q* + . Using these forms for Q± in Thtl still leaves it 
gauge invariant. Also it is clear that Thtl is real and depends only on the spatial components 
of the vector potential. Hence this is an excellent candidate for I^A), which therefore, 
following Nair, we take it to be 

^r HTL (A) 

The scheme proceeds as in the scalar theory, except that I^A) gives rise not only to a 
mass term for the free propagator, but also to higher-point interaction vertices. At one loop 
only the three- and four- point vertices are needed, and to this order the subtracting term 
uses only the quadratic contribution. Thus the gap equation reads, pictorially 




The first three graphs are as in ordinary Yang-Mills theory, with conventional vertices, 
but massive gauge field propagator (solid line); 

the first graph depicting the gauge compensating "ghost" contribution, has massless ghost 
propagators (dotted line) and vertices determined by the quantization gauge, conveniently 
chosen, consistent with the above form for the propagator, to be 

£ gaugc = ±V • A(l - /i 2 /V 2 )V • A 

fixing « \ i j / 

The remaining three graphs arise from Nair's form for the hard thermal loop- inspired I^A), 
with solid circles denoting the new non-local vertices. As it happens, the last graph with 
the four-point vertex vanishes, while the three-point vertex reads 

%t(p,q,r) = 

ifi 2 f abc ( 1 /p • q r ■ q\ r 'P/ x] 

H IT PiPjPk ~ -^-{qiqjPk + qiPjQk + Piq-flk) \ + 5 permutations 



3!(p x q) 2 [3 V p 2 r 2 jr^jr, 3r2 

p + q + r = 

The permutations ensure that the vertex is symmetric under the exchange of any pair of 
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index sets (a i p), (b j q),(c k r). Inverse powers of momenta signal the non-locality of the 
vertex. [We discuss the SU (N) theory, with structure constants f abc .} 
The result of the computation is 



Nyrab 

ij 



u YM + u N 



Hjj M is the contribution from the first three Yang-Mills graphs and IT^ sums the graphs 



from I/j,{A). The reported results ar 



Ul M (p) = N(8 ij -p i p j ) 



ij 



+ Npipj 



Npipj 



QAttll lQix I p 2/i 



\32tcll 8tc J p 2/i 8ir 



16n 
H p 



64 



32 



f3p^ + 3 E \2 ±tan _ 1 p_ 

\ 6ATCLL 167T J p 2/i 



p 



ll 



iTTLl \p z 



+ 1 



( + JL\^L tan" 1 + 

\32ttli 8tt y p 2/i 8vr 32 



16tt 8 



The Yang-Mills contribution IK is not separately gauge-invariant (transverse) owing to 
the massive gauge propagators. [At ll = 0, TiL M reduces to the standard resultllS: N(5ij — 
Pipj) (— c^p)-] The longitudinal terms in IT™ are canceled by those in IL^, so that the total 



is transverse. 



Il»(p) = N(6 ij -p i p j ] 



(-hp 2 1 \ 2/i t p 
— h —fi ] — tan — 



327T/i 27T 



V 



2/i 87T/i V p 1 J p LL &TTp Z 



[Dimensional regularization is used to avoid divergences.] 

Before proceeding, let us note the analytic structures in the above expressions, which 
are presented for Euclidean momenta, but for the gap equation have to be evaluated at the 
Minkowski value p 2 = —ll 2 < 0. Analytic continuation for the inverse tangent is provided by 



1 

— tan 

x 



1 



x 



2V^ 



In 



x/ 



ar 



Evidently IEj (p) possesses threshold singularities, at various values of —p . 

There is a singularity at p 2 = —All 2 (from tan -1 arising because the graphs in the 
Figure containing massive propagators describe the exchange of two massive gauge "parti- 
cles". Moreover, there is singularity at p 2 = —ll 2 (from tan -1 -) and also, separately in II 



YM 

ij 



and 11^, at p = (from the ±^ terms). These are understood in the following way. 
Even though the propagators are massive, the non-local three-point function contains \, 



17 



Jf, contributions, which act like massless propagators. Thus the threshold at p 2 = — /i 2 
arises from the exchange of a massive line (propagator) together with a massless line (from 
the vertex). Similarly the threshold at p 2 = arises from the massless lines in the vertex 
(and also from massless ghost exchange). The expressions acquire an imaginary part when 

the largest threshold, p 2 = 0, is crossed: Ufj M and Ilfj are complex for p 2 < 0. 

In the complete answer, the p 2 = thresholds cancel, and the singularity at the p 2 = 



/j, 2 threshold is extinguished by the factor (^ + l) 2 . Consequently IT^ becomes complex 
only for p 2 < —fi 2 , and is real, finite at p 2 = — fi 2 . 



nS(p) 



Nu 

~ PiPj) (21 In 3 - 4) 



From the gap equation in the last Figure, the result for the mass is0 

N N 

u= (211n3-4) ~ 2.384 — 

327T An 

[in units of the coupling constant g%, (or g 2 T), which has been scaled to unity]. 

Before accepting this plausible answer for /i, it is desirable to assess higher order cor- 
rections, for example two-loop contributions. Unfortunately, an estimateE-0 indicates that 79 
graphs have to be evaluated, and the task is formidable. 

An alternative test for the reliability of the above approach and for assessing the 
stability of the result against corrections has been proposed.EJ 

It is suggested that the gap equation be derived with a gauge invariant completion 
different from Nair's. Rather than taking inspiration from hard thermal loops (which after 

all have no intrinsic relevance to the three-dimensional gauge theory^ , the following formula 
for 1^ is taken 



I,(A) = /i 2 J d 3 x tr F^F* 



where D 2 is the gauge covariant Laplacian. While ultimately there is no a priori way to 
select one gauge-invariant completion over another, we remark that expressions like the above 
appear in two-dimensional gauge theories (Polyakov gravity action, Schwinger model) and 
are responsible for mass generation. If two- and higher- loop effects are indeed ignorable, this 
alternative gauge invariant completion, which corresponds to an alternative resummation, 
should produce an answer close to the previously obtained one. 

With the alternative 1^, the graphs are as before, where the propagator is still given by 
the previous expression. However, the three- and four- point vertices in I^A) are different. 
One now finds for the non-local three-point vertex 

• 2 

Ijjj Pk 

V£k(p, q, r) = -—f abC (Sijq ■ r + qffij) -r-^ + 5 permutations 

3! p z q z 



t Recall that the hot thermal loop generating functional is related to the Chern-Simons eikonal. Since 
the Chern-Simons term is a three-dimensional structure, this fact may provide a basis for establishing the 
relevance of the hard thermal loop generating functional to three-dimensional Yang Mills theory. The point 
is under investigation by D. Karabali and V. P. Nair. 
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p + q + r = 
and the non-local four-point vertex reads 



\/abcd( \ _ A* 2 f abe rede J U ^0 

v^yfei IP; q> r ) s j - ^ 1 2 J n "p2" ^2 

Pm i \ Pm ^0 , 



-e 




i£jWm - £imn£k£n-Y(P ~T - s)j + t imn e lon -^ — {j) - r - s)j(p + q - s) 



4 p 2 p^ s 2 7 ^ 

p+g+r+s=0 

These vertices do not affect the first three graphs in the Figure so that n^ M is as 
before. However, in the last three graphs the alternative non-local vertices produce the 
following result, with the help of dimensional regularization, 

% M = N( Sij - ftft) ((^ + + Jg- + 

_ f , p 4 _ p 2 , /V + x \ 2 a tan -i £ 

Y327T/X 5 167r/i 3 167r/i 327ry \p 2 J p \x 
p 2 3/i 49yU 3 /i 5 p 5 p 3 p ^ 




32?r// 16tt 96vrp 2 327rp 4 128/i 4 32/i 2 16/ 

-^M ( ( — + ~) ~ ta^ 1 ^ + ^ - ^] 
^ J \\327r// 8tt y p 2{jt 8tt 32/ 

A check on this very lengthy calculation is that summing it with Yang-Mills contributions 
yields a transverse expression. 



Uijip) = N(5ij -pipj 




Another check on the powers of - is that the above reduces to the Yang-Mills result at 
// 0.0 

Just as Nair's expression, the present formula exhibits threshold singularities: at — p 2 = 
4/z 2 , which are beyond our desired evaluation point — p 2 = /i 2 ; there are also threshold 
singularities at — p 2 = fi 2 , which are extinguished by the factor + l) 2 ; however, those at 

p 2 = do not cancel, in contrast to the previous case — indeed n^ fjo) diverges at p 2 = 0, 
and is complex for p 2 < 0. [It is interesting to remark that the last graph in the last 
Figure, involving the new four-point vertex, which vanishes in Nair's evaluation, here gives a 
transverse result with unextinguished threshold singularities at — p 2 = fi 2 and at p 2 = 0. The 
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protective factor of (p- + l) 2 arises when the remaining two graphs are added to form II 

and these also contain non-canceling p 2 
contribution.! 



threshold singularities, as does the Yang-Mills 



Although Uijip) 



is finite, it is complex and the gap equation has no solution 



for real fi 2 , owing to unprotected threshold singularities at p 2 = 0, which lead to a complex 
Uij (p) for p 2 < 0. 



a = (29±ln3 -22±) ± iN 

p 32tt V 4 3 ' 128 

N 

a ~ 1.769— 

47T 

It may be that the hot thermal loop-inspired completion for the mass term is uniquely 
privileged in avoiding complex values for — fi 2 < p 2 < 0, but we see no reason for this.f] 
Absent any argument for the disappearance of the threshold at p 2 = 0, and reality in the 
region — /i 2 < p 2 < 0, we should expect that also the hot thermal loop- inspired calculation 
will exhibit such behavior beyond the 1-loop order.^] 

Thus until the status of threshold singularities is clarified, the self-consistent gap equa- 
tion for a magnetic mass provides inconclusive evidence for magnetic mass generation. More- 
over, if there exist gauge invariant completions for the mass term, other than the hard thermal 
loop-inspired one, that lead to real U.^ at p 2 = —fi 2 , it is unlikely that they all would give 
the same fi at one loop level, which is further reason why higher orders must be assessed. 



■t-We note that Nair's hot thermal loop- inspired vertex N Vijk abc is less singular than the alternative V°j^, c , 
when any of the momentum arguments vanish. Correspondingly Hfj(p) is finite at p 2 = 0, in contrast 
to Uij(p) which diverges at However, we do not recognize that this variety of singularities at p 2 = 
influences reality at p 2 = —fi 2 ; indeed the individual graphs contributing to Ilfj are complex at that point, 
owing to non-divergent threshold singularities at p 2 = that cancel in the sum. 

§V.P. Nair states that at the two loop level, there is evidence for ln(f + ^r) terms, but it is not known 

2 

whether they acquire a protective factor of + 1). 



20 



References 



1. For a summary, see R. Jackiw and G. Amelino-Camelia, "Field Theoretical Background 
for Thermal Physics" in Banff CAP Workshop on Thermal Field Theory, F. Khanna, 
R. Kobes, G. Kunstatter, and H. Umezawa eds. (World Scientific, Singapore, 1994). 

2. R. Jackiw and S. Templeton, "How Super-Renormalizable Interactions Cure their In- 
frared Divergences," Phys. Rev. D 23 (1981) 2291. 

3. R. Pisarski, "How to Compute Scattering Amplitudes in Hot Gauge Theories," Physica 
A158 (1989) 246. 

4. E. Braaten and R. Pisarski, "Soft Amplitudes in Hot Gauge Theory: a General Analysis," 
Nucl. Phys. B337 (1990) 569; J. Frenkel and J. C. Taylor, "High- Temperature Limit of 
Thermal QCD," Nucl. Phys. B334 (1990) 199. 

5. R. Jackiw and V. P. Nair, "High Temperature Response Function and the Non-Abelian 
Kubo Formula," Phys. Rev. D 48 (1993) 4991. 

6. J. C. Taylor and S. M. Wong, "The Effective Action of Hard Thermal Loops in QCD," 
Nucl. Phys. B346 (1990) 115. 

7. R. Jackiw, "Gauge Theories in Three Dimensions at Finite Temperature)" in Gauge 
Theories of the Eighties, R. Ratio and J. Lindfors, eds. Lecture Notes in Physics 181 
(1983) 157 (Springer, Berlin). 

8. R. Efraty and V. P. Nair, "Action for the Hot Gluon Plasma Based on the Chern-Simons 
Theory," Phys. Rev. Lett. 68 (1992) 2891, "Chern-Simons Theory and the Quark-Gluon 
Plasma," Phys. Rev. D 47 (1993) 5601. 

9. D. Gonzales and A. Redlich, "A Gauge Invariant Action for 2+1 Dimensional Topological 
Yang-Mills Theory," Ann. Phys. (NY) 169 (1986) 104; G. Dunne, R. Jackiw and C. 
Trugenberger, "Chern-Simons Theory in the Schrodinger Representation," Ann. Phys. 
(NY) 149 (1989) 197. 

10. J. -P. Blaizot and E. Iancu, "Kinetic Equations for Long- Wavelength Excitation of Quark- 
Gluon Plasma," Phys. Rev. Lett. 70 (1993) 3376, "Soft Collective Excitations in Hot 
Gauge Theories," Nucl. Phys. B 417 (1994) 608. 

11. R. Jackiw, Q. Liu, and C. Lucchesi, "Hard Thermal Loops, Static Response and Com- 
posite Effective Action," Phys. Rev. D 49 (1994) 6787. 

12. P. Kelly, Q. Liu, C. Lucchesi, and C. Manuel, "Deriving the Hard Thermal Loops of 
QCD from Classical Transport Theory," Phys. Rev. Lett. 72 (1994) 3461, "Classical 
Transport Theory and Hard Thermal Loops in the Quark-Gluon Plasma," Phys. Rev. 
D 50 (1994) 4209. 

13. E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics (Pergamon, Oxford UK, 1981). 

14. J. -P. Blaizot and E. Iancu, "Non-Abelian Excitation of the Quark-Gluon Plasma," Phys. 
Rev. Lett. 72 (1994) 3317, "Non-Abelian Plane Waves in the Quark-Gluon Plasma," 
Phys. Lett. B326 (1994) 138. 

15. V.P., Nair, "Hard Thermal loops on a Moving Plasma and a Magnetic Mass Term," 
Phys. Lett. B352 (1995) 117. 

16. For a summary, see V.P. Nair, "Chern-Simons and WZNW theories and the quark-gluon 
plasma," Lectures at Mt. Sorak symposium, Korea, June 1994. 

17. G. Alexanian and V.P. Nair, "A self-consistent inclusion of magnetic screening for the 
quark-gluon plasma," Phys. Lett. B352 (1995) 435. 

18. S. Deser, R. Jackiw and S. Templeton, "Topologically Massive Gauge Theories," 
Ann. Phys. (NY) 140 (1982) 372 (E) 185 (1988) 406. 



21 



R. Jackiw and S.-Y. Pi, "Threshhold Singularities and the Magnetic Mass in Hot QCD" 
Phys. Lett. B368, 131 (1996). 



22 



